Silicon quantum computation based on magnetic dipolar coupling 
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A dipolar gate alternative to the exchange gate based Kane quantum computer is proposed where 
the qubits are electron spins of shallow group V donors in silicon. Residual exchange coupling is 
treated as gate error amenable to quantum error correction, removing the stringent requirements 
on donor positioning characteristic of all silicon exchange-based implementations [B. Koiller et ai, 
Phys. Rev. Lett. 88, 027903 (2002)]. Contrary to common speculation, such a scheme is scalable 
with no overhead in gating time even though it is based on long-range dipolar inter-qubit coupling. 
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I. INTRODUCTION 

Since the seminal exchange gate proposal of Loss and 
DiVincenzo research on semiconductor spin quantum 
computation has focused on implementations based on 
the electron exchange interaction 1 ". For silicon donor im- 
purities the use of exchange coupling is problematic since 
the exchange energy depends sensitively on donor po- 
sitioning due to the quantum interference arising from 
the six-fold degeneracy of Si conduction banc£. This re- 
sults in the necessity of donor positioning within one sil- 
icon bond (2.4 A) otherwise severe tuning requirements 
will adversely affect the scalability of this implemen- 
tation (in addition to many donor pairs having nearly 
zero exchange). This problem is attracting considerable 
attention 3 , since Si spin quantum computer architecture 
is an active research area, and donor spins in nuclear- 
spin- free silicon (pure 28 Si) are expected to have very 
long coherence timea 4 i 5 i 6 (T 2 ~ T\ ~ 10 3 s 7 ). Here we 
propose the magnetic dipolar interaction rather than the 
exchange interaction between well separated donor elec- 
tron spins as a solution to this problem. The residual 
exchange interaction is treated as a source of imperfec- 
tion in the dipolar gate, whose error probability can be 
kept below 10 -4 per operation. Hence the exchange in- 
teraction can be ignored as long as error correction is 
applied, leading to no necessity of gate tuning. This is 
possible due to the long range character of dipolar cou- 
pling (proportional to 1/d 3 , with d the inter-qubit sepa- 
ration) as opposed to the short-range nature of exchange 
[J oc d 2 5 exp (— Nevertheless this long-range char- 
acter led to speculations that a dipolar quantum com- 
puter is not scalable** We point out that this is not true, 
because only up to the 4th nearest neighbor (n.n.) cou- 
plings need to be considered, and highly efficient quan- 
tum gates can be constructed using the method of Leung 
et a/^S, which we develop further to avoid correlated er- 
ror between any two qubits inside the error correction 
manifold. Similar pulsing sequences should be useful for 



a wide variety of solid state quantum computing archi- 
tectures based on long range interactions^ 1 ". The result- 
ing architecture takes advantage of g-factor manipula- 
tion and measurement at the single spin level to avoid 
the scalability problems inherent to the ensemble NMR 
proposals (such as decreasing signal to noise ratio and 
overlapping resonances 12,13 ). 



II. DIPOLAR GATE FIDELITY IN THE 
PRESENCE OF EXCHANGE INTERACTION 

The truncated magnetic coupling between two local- 
ized spins is given by 

Hiz=w x Si z +u 2 S2z-[Di2{0,d)- J(a*,d)]S lz S 2z . (1) 

Here Si are spin- 1/2 operators, which couple to ex- 
ternal magnetic fields through the Zeeman frequencies 
uji = jiBi. Eq. (JTJ is valid provided we neglect terms pro- 
portional to Si± , which amount to a correction quadratic 
in [Dii/lyix — u> 2 )^- Hence a strong inhomogeneous 
field is needed (or inhomogeneous gyromagnetic ratios 
7i), but since Z? 12/74 < 0.01 G (see below) field differ- 
ences on the 10 — 100 G range are sufficient. The dipolar 
interaction is given by 



D 12 (6,d) = ^fi (3 cos 2 9 -l) 
d J v ' 



(2) 



where d is the inter-qubit distance and is the angle 
between the external magnetic field and the line joining 
the spins. The optimum dipolar architecture assumes 9 = 
0, e.g., an array of spins directed along B (this optimal 
coupling is denoted D\ 2 below). Eq. Q has a striking 
property: If |cos#| = l/y/3, the interaction is exactly 
zero. Hence in an array of spin qubits dipolar coupling 
can be completely suppressed as long as ±B makes one of 
the "magic" angles with the array: 9 = 54.74°, 125.26°. 
Exchange based proposals 1 usually require the donors to 
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be pushed outside the array to switch on the exchange 
interaction. In addition, 2D and 3D arrangements have 
been considered, particularly to optimize quantum error 
correctioni£ii&. In this case it may be impossible to find a 
geometry where all bonds are making a magic angle with 
the B field. Then if dipolar interaction is to be ignored, 
we will loose track of the spin state within 10 /zs (this 
time should be considered short if compared with other 
mechanisms such as nuclear spectral diffusion^). Hence 
dipolar coupling may have to be taken into account even 
in exchange gate quantum computing architectures. For 
the exchange interaction between two hydrogenic donors 
we use the asymptotic expression 8 



J{a*,d) 



1.6 e 2 

he a* 



5/2 



exp 



(3) 



valid for inter-donor distances d much larger than the 
Bohr radius a*. Eq. 10) is to be regarded as an enve- 
lope for the strong oscillations of the exchange energy 
stemming from conduction band degeneracy 2 * The Bohr 
radius a* is related to the experimental donor ground 
state energy Ed, see Table U^. 

We will now show how a silicon donor quantum com- 
puter can be implemented with the dipolar interaction 
and single spin rotations. The effect of the exchange in- 
teraction will be treated as an error, leading to a lower 
bound on qubit separation. Assuming J = in Eq. l[TJl. a 
controlled-Z (CZ) gate is obtained by free evolution dur- 
ing the time interval tqz = tt/-Di2 together with g- factor 
shifts^, 



U C z = e-- l e- lbl *e-? lb2 * exp \~i—Ui 2 j . (4) 
Below we show how to correct for the Zeeman frequencies 



Donor E d [meV] a* [A] d opt [A] d [A] t C z [fis] 



Sb 


43 


18.6 


315 


263 


150 


P 


45 


18.2 


307 


256 


140 


As 


54 


16.6 


279 


232 


105 


Bi 


71 


14.5 


241 


200 


68 



TABLE I: A group V donor electron spin quantum computer, 
where free evolution of the spin-spin dipolar interaction im- 
plements CZ gates. Here we show donor electron ground state 
energies Ej~, Bohr radius a*—-, optimum inter-qubit dis- 
tance d pt (for the exchange interaction to be ignored within 
10 -4 error probability), inter-qubit distance do (such that 
D12 = J), and the CZ gate times. Fastest gate times are 
obtained for bismuth donors. 



E = 1 — F 2 ). If one wants to ignore exchange interac- 
tion, all that needs to be done is to keep a 2 /16 less than 
some critical bound p, for example the 7-qubit encoding 
threshold p = 10~ 419 . Hence we have J ' jDyx < ^fp, or 



J_ 

D 



0.02A 



2 , d x 11/2 



exp -2— < 1CT 2 , (8) 



for p — 10 4 . The leng th scale for the prefactor in this 
expression is given by ■ s /2e/lM~//e » 0.02 A. The range 
where this inequality is satisfied is approximately given 
by < d < 0.03a* and d > 17a*. The first condition 
arises due to the divergence of the dipolar interaction, 
and is not useful here [also Eq. |(3J is only appropriate 
for d a*}. The physical solution is the second one, 
which is optimal (fastest gate) for d opt « 17a*. Tabled 
shows d op t for various donors together with their CZ gate 
time [7r/Di2, see Eq. Q]. 



We now search for the minimum inter-qubit distance d 
so that J can be ignored. A residual exchange interaction 
J will add an additional evolution operator to Eq. @ 



U(a) = exp {-iaS lz S 2z ), 



(5) 



with a — ■kJ/D\2. This causes phase error in the CZ 
gate, which is better evaluated by looking at two input 
states orthogonal to each other. Equivalently we look at 
the CNOT gate (obtained by a basis change on the CZ, 
?7cnot = e~ l ^ S2y Ucze z: % S2v )- Therefore the "erroneous" 
evolution is given by U(a)UcNOT, where 



U(a) 



U{a)e i 



(6) 



is a 4 x 4 matrix with elements equal to cos (a/4), 
sin (a/4), and 0. The error due to a finite a can be 
evaluated by calculating the fidelity functions 



F{\i>),a} 



I ^CNOT U(a) f/cNOT 1 1p) 



(J) 



which are simply given by | cos j \ , leading to an error 
probability of a 2 /16 for small a (the error is given by 



III. A SCHEME FOR DECOUPLING 
LONG-RANGE INTERACTIONS 

Up to now we have shown that dipolar coupling be- 
tween two donors can generate precise two-qubit evo- 
lution i.e. a dipolar coupled-qubit Si gate can be con- 
structed. However the situation becomes complicated 
when we consider an array of many donors. Particularly 
the long range nature of the dipolar interaction implies 
every spin in the array will be coupled to each other, 
raising questions about the scalability of this proposal 
(This was one of the original motivations for introducing 
the exchange gate since exchange can be exponentially 
suppressed by electrically controlling wave function over- 
lap). For example, it is possible that the complexity of 
the pulsing sequences (leading to the desired quantum 
algorithm) might scale exponentially with the number of 
qubits, effectively making the problem of determining the 
evolution as hard as any mathematical problem a quan- 
tum algorithm is constructed to solve 9 . Nevertheless this 
is not true for the case considered here, because using the 
same argument leading to the discard of the exchange 
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interaction we can neglect (within the 10~ 4 threshold) 
dipolar coupling between any spin and its 5th or higher 
n.n. [By Eq. JBJ, D lk /D 12 = l/k 3 , which is less than 
1CP 2 for k > 5]. Hence Eq. generalized to a ID spin 
array is 



repetition), 



H 



i+4 

£ 



Dij Si z S j z j 



(9) 



where i is an integer labeling the location of each donor 
(i is assumed positive as well as negative). The finite 
coupling range allows us to develop quantum gates using 
a sequence of 7r-pulses applied to subsets of the spins 
[each 7r-pulse is given by X^> — exp(inSi x )]. The key 
point is that the interaction between any two spins can 
be canceled using two 7r-pulsesi^, 



exp 



-i-DSi z S 2 



; )*( 2 )exp(- 



-t-DS lz S 2z )X^ =1, 
(10) 

where / is the identity operator. Our task is now to 
find the "decoupling" scheme which completely refocuses 
Eq. © after some time interval r (therefore enabling 
single qubit rotation on any spin) and also to produce 
sequences for "selective recoupling", which provide CZ 
evolution for any n.n. pair. For this task we use the 
method of Ref. [jiJ which consists in constructing sign 
matrices S n representing the 7r-pulses. A n x m sign 
matrix has each element equal to ±1 (denoted simply by 
±), and correspond to a system of n spins where evolution 
during a time r is divided into m time intervals. If spin 
i has its interaction reversed in any Z-th time interval [by 
application of X^> before and after this time interval], 
then (S n )u = — 1, otherwise (S n )u — +1. For example, 
Eq. 1|10|) corresponds to 



(11) 



The interaction between two spins i, j is decoupled if the 
rows i and j of S n disagree in sign for half of the m 
time intervals. Equivalently, the inner product between 
these rows is zero. This property leads to a connection 
with the theory of Hadamard matrices: H w isanxn 
Hadamard matrix if and only if ■ (Hn) T = nl . Hence 
a possible solution for the decoupling problem of n spins 
is to construct S n from n rows out of a Hadamard Hn 
where n > n. Actually such a solution turns out to be 
the most efficient one (the smallest n satisfying n > n is 
the minimum number of intervals m in the set of possi- 
ble S n ) because one can show that it is impossible to add 
an additional row orthogonal to . Hadamard matrices 
exist for n = 1, 2 [Eq. (H])]. 4, 8, 12, . . . (see Ref. and 
references therein). The finite coupling range of Eq. © 
suggests n — 12 as a convenient solution. S n can be as- 
sembled asanx 12 matrix composed of 7 ordered rows 
from H12 (identical rows are 7 rows apart). For the par- 
ticular case of n = 14 (general n is obtained by row 



Si 




(12) 



Here we extract the first (and last) 7 rows of Su from 
H12 with the first row (+ + ...) excluded so that Zee- 
man splitting is also canceled. £14 requires a total of 80 
7r-pulses which are applied in 12 sets (less than 14 pulses 
are applied in each set because X^ 2 = I - hence no ro- 
tations need to be applied when the sign is the same for 
neighboring time intervals). An array of n spins will re- 
quire less than 6n pulses. Selective recoupling is achieved 
by choosing identical rows for the spins which are to be 
coupled. These rows are chosen from the 4 remaining 
rows of -Hi 2, for example 



/ + + + + 



o 1 1 



+ 
+ 



+ 
+ 



+ - + - 
+ - + - 



- + + 

- + + 



+ 
+ 



v + 



+ 
+ 



- + - + - + - + 

- + - + - + - + 



(13) 

implements CZ operations between spins 6, 7 and 10, 
11 in parallel (bold). We point out that each 7-qubit 
structure in Eqs. i|12|) and l|13|l form an error correction 
block for the Steane codeA2>. Note that residual dipolar 
interaction couples qubits in different blocks. This is im- 
portant for the assumption of uncorrelated errors within 
each block and the validity of the 10~ 4 threshold (see 
section IV) . The spurious couplings lead to error of the 
order of 7~ 6 - 10~ 5 in Eq. |Q2J and 22~ 6 - 10~ 8 for 
selective recoupling when all blocks execute CZ in par- 
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allel (massive parallelization of CZ is needed for efficient 
computation and quantum error correction). 

Therefore the complete gate time for a large ID array 
is the same as for two donors (approximately 100 fj,s, be- 
ing optimal for bismuth - see Tabled . This shows that a 
dipolar donor electron spin quantum computer is reliable: 
If the silicon lattice is isotopically purified (free of 29 Si 
nuclear spins), the coherence time will be limited by the 
spurious exchange and dipolar couplings, with a quality 
factor of the order of 10 4 . A key advantage of this archi- 
tecture is the inter-qubit distance, which is three times 
larger than other proposals for donors 1 . Also there is 
no need for an inter-qubit "J" gate, or any electrical con- 
trol over wave function overlap 1 ^. This should make gate 
lithography much simpler (one needs to incorporate g fac- 
tor control 21 and single spin measurement/initialization 
electrodes 2 ^ on top of each donor). 

The considerations above can be generalized to any 
long-range coupling D oc l/d r . The number of n.n. 
which need to be decoupled is given by Max (k) < p~ l / 2r , 
where p is the desired error probability. Hence rcz needs 
to be broken into n ~ p^ 1 / 21 " time intervals. For example, 
r > 1 and p = 10~ 4 leads to n ~ 100. Implementation 
of any quantum gate is possible as long as the time for 
single spin rotation is much less than Tcz/n. The dipolar 
case considered here clearly satisfies this criterion, since 
i~cz/n ~ 10 £*s (rotation times of the order of 0.1 [is are 
easily achievable). Finally, notice that this approach for 
decoupling can also be applied to general anisotropic ex- 
change interactions, since these can be transformed into 
the SizSj Z form by appropriate spin rotations. 



IV. ERROR CORRECTION OF RESIDUAL 
LONG-RANGE COUPLING 

Here we show how imperfections arising from spuri- 
ous long range couplings connecting qubits in distinct 
quantum error correction blocks can be corrected by the 
usual syndrome diagnosis (projective measurement on 
each block The proof presented here is based on the 
simplest error correction code, the "three bit flip code" 
(section 10.1.1 of Ref . 123^ . However, we emphasize that 
these results are easily extended to the complete 7-bit 
Steane's code 1 ^ which correct for any type of continuous 
error on each qubit within its block. The essence of our 
proof is that the syndrome measurements on each block 
effectively destroys error correlation between qubits be- 
longing to different blocks. 

Consider two error correction blocks constituted by 
qubits 1,2,3 (first block) and 4,5,6 (second block). The 
residual coupling Hamiltonian is 

H = -Ac 1 {S lz S iz + S 2z Sr oz + S 3z S 6z ) , (14) 



and the evolution operator after one "clock time" r is 

U (r) = exp (—iTLr) = cos 3 (c)I + i sin (c) cos (c) 

x [o\ z o\ z + <r 2z o- 5z + o- 3z a 6z ] - cos (c) sin 2 (c) 
x [o-i z a 2z o- 4:Z a 5z + o-x z a 3z a^ z G^ z + o-2 Z a 3z a 5z a 6z ] 
-i sin 3 (c) [a lz a2 Z o- 3z a4, z o- 5z a 6z } . (15) 

Here c = c'r = ttDu/D^ is much less than one (in the 
case of 7-qubit blocks c ~ 7 -3 ). To map this problem into 
the bit-flip code we use the y basis for our spin qubits: 

10) = l + y> = ^=(lT>+i|i», (16) 
|i> = I -y) = ^ (IT)-* ID). (17) 

In this basis the Pauli matrices of Eq. I|15fl act as a bit 
flip operator (a z \Q) = |1), tr a |l) = |0)). Eq. ((T5J con- 
tains three contributions: (1) simultaneous bit flips of 
one spin in block 1 and another in block 2. This leads 
to error probability of ~ c 2 for each spin in each block 
[square of the amplitude, see Eq. Q]. Note that error is 
correlated between blocks; (2) Simultaneous double bit 
flip in both blocks, with probability equal to the square 
of the single bit flip probability (~ c 4 ). This process has 
the same order of magnitude of two independent single 
bit flip errors occurring at the same time. This type of 
error is only corrected after two concatenations of error 
correction are in place; (3) three bit flip error in both 
blocks, probability is cube of single bit flip, equivalent to 
three simultaneous independent bit flips. 

We start by assuming blocks 1 and 2 store the state 

|^(0)) = oi|00) + <z 2 |01) + a 3 |10) + a 4 |H). (18) 

For fault tolerant quantum computing, we add 4 addi- 
tional ancilla qubits which encode the state as 

|V>(0)) = ai|000,000) + a 2 |000, 111) 

+a 3 |lll,000) + a 4 |lll,lll). (19) 

Time evolution under the spurious coupling Hamiltonian 
[Eq. dSJ] yields 

|*(r)) = cos 3 (c)|*(0)) - isin 3 (c)\<p) + i sin (c) cos 2 (c) 
x {ai [| 100, 100) + i tan (c) 1 01 1,011)] 
+a 2 [|100,011) + £tan(c)|011,100>] 
+a 3 [|011,100) + £tan (c)|100, 011)] 
+a 4 [|011,011) + itan(c)|100,100)]} 
+ •••, (20) 

where the swapped state \<f>) is given by 

\(t>) = oi|lll,lll) + a 2 |lll,0Q0) 

+a 3 |000,lll) +a 4 |000,000). (21) 
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Error correction proceeds with projection measurements 
over the syndromes 0,1,2, and 3 in each block k = 1, 2: 

F (fc) = |000)(000| + |111)(111|, (22) 

p[ k) = |100)(100| + |011)(011|, (23) 

F 2 (fc) = |010)(010| + 1 101) (101| , (24) 

P 3 (fc) = |001)(001| + 1 110) (110|. (25) 

Depending on the outcome of the measurement we apply 
the corresponding correction operator u[ k ^ (for example, 

t/« = J, U[ 1] = a u , u£\ = a 2z , =a 3z ). The final 
corrected density matrix is an incoherent superposition 
of each possible error: 

id 

x(*(t)\P^U^P^U^. (26) 

Because error is correlated between blocks the projec- 
tion P^P^ onto state |^(r)) is zero unless i — j. For 
example, 

2 

n^o (fc) ^o fc) l*W> = cos 3 (c)|*(O))- lS in 3 ( C )|0), 
fe=i 

2 

TJ[|WpW|$( T )) = isin(c)cos 2 (c) 
fc=i 

x[|tf(0))+itan(c)|$], (27) 

for / = 1,2,3 [Note that Eq. (JJOJ omitted the syndrome 
subspaces 2, 3]. Finally it is straightforward to calculate 
the fidelity squared, 

|F| 2 = (tt(0)|p c ||*(0)) 

= cos 4 (c) [1 + 2 sin 2 (c)] 

+ sin 4 (c) [l + 2cos 2 (c)] |(*(0)|^)| 2 . (28) 

Maximum error occurs when the second term of Eq. I|28|) 
is zero. This leads to 

Max (E) = 1- Min (|F| 2 ) » 3c 4 + C(c 6 ), (29) 

which is the square of the error without error correction. 
Hence neglecting small dipolar coupling between differ- 
ent error correction blocks is for all practical purposes 
equivalent to having a source of independent uncorre- 
cted error identical to the one assumed in the quantum 
error correction literature. 



V. DISCUSSION 



We now consider the feasibility of our dipolar QC pro- 
posal for III-V semiconductor donor impurities and quan- 
tum dots. Although these materials have a small effective 
mass (implying higher a* and <i pt)> some of the nar- 
row gap semiconductors have quite large bulk g factors, 
enhancing dipolar coupling. A simple estimation is ob- 
tained from the relation rcz ~ (0.3/ro*) 3 (2/q) 2 x 100 /is. 
Using the parameters of Ref. |24| we get tqz ~ 0.1 s for 
GaAs and tqz ~ 1 ms for GaSb, InAs, and InSb donors 
impurities (quantum dots have dipolar gate times higher 
by approximately a factor of 10 due to larger Bohr radii). 
Hence our proposal is not feasible for GaAs, but might 
work for the narrow gap III-V materials as long as deco- 
herence due to nuclear spectral diffusion is suppressed by 
nuclear polarization^. In this case spin-flip followed by 
phonon emission will be the dominant decohering pro- 
cess. Adjusting the external magnetic field, coherence 
times of the order of a few seconds are achievable 2 ^, sug- 
gesting the possibility of quality factors greater than 10 3 
in a narrow gap donor dipolar quantum computer, which 
does not require exchange interaction control and can be 
constructed with current lithography techniques. 



In conclusion we consider a quantum computer archi- 
tecture based on dipolar-coupled donors in silicon. Al- 
though gate times are considerably longer than exchange- 
based implementations (albeit same time scales as the 
solid state NMR proposals^), one does not need atomic 
precision donor implantation or electrical control of two- 
qubit couplings. Particularly "top-down" construction 
schemes based on ion implantation should benefit from 
our proposal, because these lack precision in donor posi- 
tioning in addition to creating interstitial defects 25 (dipo- 
lar coupling is nearly insensitive to electronic structure). 
Our proposal for decoupling of short-range "always on" 
interactions together with error correction of the remain- 
ing long-range couplings apply equally well to any solid 
state implementation based on other types of long-range 
interactions (as long as the coupling is bilinear 
opening the way to implementations which do not have 
severe lithography requirements. We acknowledge useful 
discussions with B. E. Kane, J. Kempe, T.D. Ladd, T. 
Schenkel, J. Vala and W. Witzel. This work is supported 
by ARDA, LPS, US-ONR, and NSF. 
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